A New Wavelet Method for Identification of Eddies and Assessment of Incidents on Islands of the Eastern Caribbean

Using two dimensional continuous wavelet transforms, a novel method for identification of mesoscale eddies is presented to facilitate extraction of characteristics for area, amplitude, type, and location from maps of sea level anomalies. In comparison with the previously established growing method for eddy identification, it is found that the wavelet method identifies more than twice the number of eddies and is particularly better at resolving small eddies down to the 0.25 degree resolution of the data. Such research into eddy identification and tracking is significant to the assessment of eddies with potential to impact on coastlines of small islands. The method is applied to the identification of eddies on tracks towards islands of the Eastern Caribbean over 23?years. Spatial and temporal variation in rate of occurrence and magnitude is established. For Barbados there is an average of 9 anticyclonic incidents a year with maximum amplitude of typically 0.22?m in the dry seasons and 0.16?m in the wet seasons. Seasonal variation is reversed for the other islands with twice the number of anticyclonic incidents having maximum amplitudes of about 0.20?m annually.     

复杂地形景区三维导览图的设计与实现

为满足复杂地形景区对三维地图导览的需求,构建更为清晰的三维地图模型,提供良好的三维导览地图设计方案显得尤为重要。为了更好地与卫星影像进行贴合,DEM数据的采样密度要与卫星影像的分辨率一致。本文通过对比当下常用的四种空间插值方法的适用范围与运用特点,设置相关阈值及权重,直观比较了四种插值方法产生的插值结果,选择更适合复杂景区DEM插值的插值方法并对其进行精度提升,以此得到符合精度要求的DEM数据。最后利用相应地区的卫星影像进行地图投影及影像贴图,两种数据结合构建复杂景区的三维地图模型,给游客带来更为精确直观的定位信息和空间要素信息。     

Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

沉积物中硬质粘土层的研究进展

大量的钻孔资料表明,平原地区及其海域沉积物中普遍有硬质粘土层的存在。硬质粘土层在工程地质、地理学研究等方面具有重要的价值,它不仅是高层建筑的持力层,而且是研究全球环境变化事件、古气候演变的良好载体。该文阐述了硬质粘土层的一般特征、成因以及在粒度、地球化学元素、年代学、沉积环境等方面的研究进展,并展望了硬质粘土层在年代、物质来源和沉积环境等方面的研究方向,认为黄河三角洲硬质粘土层研究较少,需要开展大范围的、系统深入的、多学科的宏观和微观综合分析,挖掘硬质粘土层蕴含的古环境意义和应用价值。     

杭州湾中部实测波浪特性分析

利用杭州湾中部一年实测波浪资料,分析该地区的波参数统计变化特性,采用最小二乘法拟合得出相关参数之间的关系式,并统计分析了频谱特性。结果表明:杭州湾中部以小浪和轻浪为主,在各月分布较为均匀,年平均周期为2.97 s;最大一次波浪过程为冷空气影响所致,影响强度大于台风莫拉克;常浪向分布在东北至东南向,其中又以小浪出现较多,对应的波周期以2~4 s为主;强浪向主要分布在西北和东北方向,对应的波周期主要分布在3~5 s;多数特征波参数之间相关关系较好;由谱分析得知风浪占多数,且以单峰为主。研究结果可为相关工程和理论研究提供基础资料。